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1. Introduction 



In this paper, we investigate the plane flow of ideal plastic materials [21 El S] modelled 
by the hyperbolic system of four partial differential equations (PDE) in q = 4 dependent 
variables a, 6, u, v and p = 2 independent variables x and y, 

(a) o~ x - 2k (6 X cos 26 + 6 y sin 26) = 0, 

(b) o- y - 2k (6 X sin 26 - 6 y cos 26) = 0, 

(c) (uj, + t^) sin 26 1 + (m x — v y ) cos 26 1 = 0, 

(d) u x + v y = 0, 

where a x = da/dx, etc. The expressions ((Ha), ((Ub), are the equilibrium equations 
for the plane problem. In other words they are the Cauchy differential equations of 
motion in a continuous medium where we consider that the sought quantities do not 
depend on z. These two equations involve the dependent variables o and 6 that define 
the stress tensor; a is the mean pressure and 6 is the angle relative to the x axis in 
the counterclockwise direction minus 7r/4. The equation ([Uc) corresponds, in the plane 
case, to the Saint- Venant- Von Mises plasticity theory equations, where u and v are 
respectively the velocities in the x axis and y axis directions. Moreover, we assume 
that the material is incompressible and hence that the velocity vector is divergenceless. 
This explains the presence of the equation ([Ud) in the considered system. The positive- 
definite constant k is called the yield limit and it is associated with the plastic material. 
Without loss of generality we assume that k = 1/2 (this is the same as re-scaling the 
pressure a). 



In a recent work (18], the concept of homotopy of two functions has been used to 
construct two families of exact solutions for the system formed by the two first equations 
in (CQ). In the papers [SHE], the Nada'i solution [7J for a circular cavity under normal stress 
and shear and the Prandtl solution [8j for a bloc compressed between two plates has been 
mapped by elements of a symmetry group of the system consisting of ((Ua) and ([Ub), 
in order to calculate new solutions. In addition in |9], simple and double Riemann wave 
solutions for the system ([1]) were found using the method of characteristics. However, 
as is often the case with this method, those solutions rely on numerical integration 
for obtaining the velocities u and v. Symmetries of the system ([!)) were found in [lj. 
However, the Lie algebra of symmetries, was not complete because of the absence of the 
generators B\, B3, B4, B±, Bq and K (or equivalents) defined by (jlj) in Section [2] of this 
work. Moreover, we found two infinite-dimensional subalgebras spanned by X± and X 2 
provided below in equation ([5]). The generator X\ was known [TJ as a symmetry of the 
two first equations in ([T|), but it is shown in this paper to also be a symmetry of the 
complete system (CQ). The generator X2 is a new one. A classification of one-dimensional 
subalgebras was performed in [TJ. To our knowledge, no systematic Lie group analysis 
based on a complete subalgebra classification in conjugacy classes under the action of 
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the symmetry group G of the system ([I]) that includes the new found generators has 
been done before. 

The goal of this paper is to systematically investigate the system (OQ) from the 
perspective of the Lie group of point symmetries G in order to obtain analytical solutions. 
That is, we obtain in a systematic way all invariant and partially invariant (of structure 
defect 5 = 1 in the sense defined by Ovsiannikov |10| ) solutions under the action of G 
which are non-equivalent. Invariant solutions are said to be non-equivalent if they cannot 
be obtained from one another by a transformation of G (the solutions are not in the 
same orbit). In practice, we apply a procedure developed by J. Patera et al. [TT | [T2 | [T3] 
that consists of classifying the subalgebras of £ associated with G into conjugacy classes 
under the action of G. Two subalgebras Ci C C and C\ C C are conjugate if GCiG = 
For each conjugacy class, we choose a representative subalgebra, find its invariants 
and use them to reduce the initial system (OQ) to a system in terms of the invariants 
which involve fewer variables. We illustrate these theoretical considerations with some 
classes of algebraic solutions, some of them in closed form and others defined implicitly. 
Thereafter, we draw for some solutions and specific choice of parameters the shape of the 
corresponding extrusion die. The applied method relies on the fact that the contours 
of the tools must coincide with the flow lines described by the velocities u and v of 
the solutions of the problem. For applications, it is convenient to feed material into 
the extrusion die rectilinearly at constant speed. So, the tools illustrated in this paper 
were drawn considering this kind of feeding. Based on mass conservation and on the 
incompressibility of the materials, we easily deduce that the curve defining the limit 
of the plasticity region for constant feeding speed must obey the ordinary differential 
equation (ODE) 

dy_ = V -v(x,y) _ 
dx U — u(x, y) ' 

where Uq, Vq are components of the feeding velocity of the die (or extraction velocity at 
the output of the die) respectively along the x-axis and y-axis. One should note that 
the conditions (J2J) are reduced to those required on the limits of the plasticity region in 
reference [9] when Vq = and that the curves defining the limits coincide with slip lines 
(characteristics), that is when we require 

dy/dx = tan 6(x, y) or dy/dx = — cot 8(x, y). (3) 

Thus the condition ([2]) can be viewed as a relaxation of the boundary conditions given 
in [9]. The reason we can use these relaxed conditions is that we choose the contours 
of the tool to coincide with the flow lines for a given solution rather than require the 
flow of material to be parallel to the contours. Using these relaxed conditions, we can 
choose (in some limits) the feeding speed and direction for a tool and this determines 
the limits of the plasticity region. 
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This paper is organized as follows. In section [2] we give the infinitesimal generators 
spanning the Lie algebra of symmetries C for the system ([I]) and the discrete 
transformations leaving it invariant. A brief discussion on the classification of 
subalgebras of C in conjugacy classes follows. Section [3] is concerned with symmetry 
reduction. It describes how the symmetry reduction method (SRM) has been applied 
to the system ([I]) and the method for finding partially invariant solutions (PIS). Some 
examples of solutions are presented, including invariant and partially invariant ones. We 
conclude this paper with a discussion on the obtained results and some future outlook. 

2. Symmetry algebra and classification of its subalgebras 

In this section we study the symmetries of the system ([I]). Following the standard 
algorithm j 15 j . the Lie symmetry algebra £ of the system has been determined. Using 
the notation d x = d/d x , etc., the Lie algebra of symmetry is spanned by the fourteen 
infinitesimal generators 

-Pi — d x , P 2 = d y , P3 = d u , P4 = d v , P5 = d a , 
L>i = xd x + yd y + ud u + vd v , D 2 = xd x + yd y - ud u - vd v , 
L = —yd x + xd y — vd u + ud v + dg 
B 1 = -vd x +ud y , B 2 = yd u -xd v , 
B 3 = (tr + ±sin26)d x - ±cos2# d y , B 4 = -\ cos20 d x + {a - (1/2) sin20)^, (4) 
B5 = (cr — I sin 28)d u + | cos 26 d v , B & = ^ cos 26 d u + (<t+ § sin 26)d v , 
K =(\xco$26 -y(a + \sm26))d x + {{a - \ sm26)x + \ycos26)d y 

v(|sin26 l — cr))d u + ((a + |sin2#)u — |v cos 26)d v 

+ 6d a + ad e , 
and the two generators 

X 1 = 6)d x + ^(a, 6)d y , X 2 = 0(a, ^)9 U + Vfa 0)a„, (5) 
where the coefficients ^ and must satisfy the two quasilinear PDEs of the first order 
= cos 26£ e + sin 26r] B , ^ = cos 26£ a + sin 20r/ CT , (6) 

while the coefficients <fi and ip must solve the two PDEs, of the same type as the previous 
ones, 

(j)a = - (cos 264> e + sin 26ip e ) , = - (cos 26>0 CT + sin 26i/; a ) . (7) 

Note that Xi and X 2 span infinite subalgebras. The generators D\ and D 2 generate 
dilations in the space of the independent variables {x, y} and dependent variables {u, v}. 
The generator L corresponds to a rotation. Moreover, Bi, i = 1, ... ,4 are associated 
with a type of boost and the P i: i = 1, . . . , 5 generate translations. Of these fourteen 
generators eight were already known but six are original to this paper. The newly 
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found generators are B\, B 3 — Bq and K. In addition, the system ([I]) admits two infinite 
dimensional subalgebras. The one spanned by X± was known [1] as a symmetry of the 
two first equations of the system (0Q), but it still a symmetry of the entire system (CQ). 
The infinite subalgebra corresponding to X 2 is original to this paper. The commutation 
relations for the 14— dimensional Lie subalgebra 

£ = {Pi, D 2 , B 2 , L, D u B 3 , P 4 , B 5 , B 6 , P x , P 2 , P 3 , A, P 5 } , (8) 

excluding K, are shown in table [TJ Note that the generators Pi, P 2 , P 3 , P 4 are just 



Table 1. Commutation relations for the algebra C excluding K in equation (TJ|. 
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Bi 


D 2 


B 2 


Dx 
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B 3 


Bi 


B 5 


Be 
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p 3 


Pi 
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-D 2 
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particular cases of X\ in ([6]) while P3, P4, P5, Bq are particular cases of X 2 . Nevertheless, 
we include them in the classification under conjugacy classes that we consider in the 
sequel. We exclude K from C because it cannot span a finite Lie algebra together with 
Pi, P 2 , P3, P4, B 3 , P 4 , B 5 , B 6 . The maximal finite symmetry Lie algebra that includes K 
is 

S = {B 1 ,D 2 ,B 2 ,K,L,P 5 ,D l }. (9) 
The commutation relations for S are given in table [2 One should note that the system 



Table 2. Commutation relations for the algebra S. 
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Bi 


D 2 


B 2 


K 


L 


p 5 




Bi 





2Bi 


-D 2 
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D 2 







2B 2 














B 2 


D 2 


-2B 2 
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-p 5 


-L 





L 











P 5 











P 5 











L 











Di 
























([T]) is also invariant under the discrete transformations: 

Pi : x 1— > —x, y I—)- —y, a a, 9^9, «4 u, v M- v; 
R 2 : x i-> x, y t-ty, o >-)■ a, 9^9, ti4 —u, v >-)■ —v. 
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These transformations Ri and R 2 are rotations by an angle it in the plane of independent 
variables x, y and of dependent variables u, v respectively that induce the automorphisms 
of the Lie algebra £: 

K X :D X ^ D h D 2 ^ D 2 , Bi ^ -Bi, 1 = 1..4, B t i-> B u i = 5, 6, 
Pi^ -P u i = l,2, Pi ^ P h i = 3,4,5, Lh>I, K^K; 

K 2 :D X ^D U D 2 ^ D 2 , B, t ^ -B h i = 1,2, 5, 6, B { h-> B h i = 3, 4, 
P ^ P, z = 1,2,5 Pi^-Pi, i = 3,4, L^L, K ^ K. 

Since we look for solutions that are invariant and partially invariant of defect structure 
5 = 1, we only have to classify the subalgebras of codimension 1 and 2. We consider 
separately the two distinct finite dimensional symmetry Lie algebras discussed above. 

2.1. Classification into conjugacy classes of the subalgebra L. 

In order to apply the method of classification into conjugacy classes developed in [TT|[T2]. 
one first chooses a decomposition of the Lie algebra. We now describe the decomposition 
used for the classification of algebra C. We begin by factoring C into the semi-direct 
sum of the one-dimensional subalgebra {P5} and the ideal 

M = {Pi, D 2 , B 2 , L, D u P 3 , P 4 , P 5 , B 6 , P u P 2 , P 3 , P 4 } , 

i.e. 

£ = {P 5 }>M. (12) 

Before we can apply the classification procedure to the semi-direct sum ( fl~2]) . we have 
to know the classification of the subalgebras {P5} and M. The subalgebra {P5} is a 
representative of its conjugacy class. In order to classify M, consider the following 
decomposition: 

M=T>M, (13) 

where 7 = {B X ,D 2 ,B 2 ,L,D X } and M = {P 3 , P 4 , P 5 , P 6 , Pi, P 2 , P 3 , P 4 } is an Abelian 
subalgebra. The subalgebra F is further decomposed into the direct sum 

F = A®11, (14) 

of a simple algebra A = {Bi,D 2 ,B 2 } and an Abelian algebra 1Z = {L,Di}. The 
classification of the simple algebra was carried out by Patera and Winternitz in their 
work [16] on the classification of subalgebras of real Lie algebras in three and four 
dimensions. For the classification of subalgebra R, the conjugacy classes are simply 
its subspaces since R is Abelian. We then use the Goursat twist to obtain a list of 
representative subalgebras for the conjugacy classes of F. The subalgebra iV is Abelian, 
so, once again, all subspaces are representative subalgebras. Therefore, any further 
decomposition may seem superflous. However, since we factored out the infinitesimal 
generator P 5 and it does not appear in M, subalgebras of type B © {0} and {0} © V, 
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where B = {B 3 , B A , B 5 , B G } and V = {Pi, P2, -P3, -P4}, have the same algebraic structure 
in M, see table [TJ That is, they have the same commutation relations with elements of 
F. Therefore, taking the decomposition 

M = B © V, 

a classification of the splitting subalgebras of the form Ti > {B © {0}}, where J 7 , C J 7 
is a representative subalgebra of the classification of J 7 , will lead to that of the splitting 
subalgebras of form J*j D> {{0} © V} when (B 3 , _B 4 , B 5 , B 6 ) is relabelled (Pi, P 2 , P3, Pa). 
Since a splitting subalgebra is a subspace of form JFi x Bi in which Bi C B is an ideal, the 
splitting subalgebras with non-zero component in B and in P are classified as follows. 
For each splitting subalgebra of form Ti > © {0}} whose basis vectors Bj are labeled 
b s , s — 1, . . . ,n — dim Hj , one adds to each basis vector b s an element p s of V under the 
most general form p s = Ylt=i VstPt, A*st £ Next, we require that the space generated 
by {b s + p s }, s = 1, . . . , n be an ideal in Pj x {6 S + p s }. This leads to constraints on the 
parameters /j, st . Once these conditions have been satisfied, one reduces the range of these 
parameters as much as possible by conjugation with Nor(Fi > {Bj © {0}} ; exp Ai). This 
completes the classification of the splitting subalgebras of M. From the list of splitting 
subalgebras, we use the procedure described in [17] to find the subalgebras which are not 
conjugate to any splitting subalgebras, the so-called nonsplitting subalgebras. Finally, 
we have to classify the semi-direct sum {-P5} >A4, which we do by using the semi-direct 
sum as developed in (TTJ [12]. Next, we reduce the range of the parameters as much as 
possible using the algebra automorphisms f TTTT) . The result of the classification of the one- 
dimensional subalgebras is presented in table [3j The classification of conjugacy classes 
for two-dimensional subalgebras consists of a list of 250 representative subalgebras. This 
list can be found in the Appendix. 



£1,1 = 


= {Bi\ 




£1,19 


= {D 2 - 


-Pi 4 


B 3 + 


aP 2 +SP 5 \ 


£1,37 


= {P 3 4 


P54- 


ePi + aP 2 } 


£1,2 = 


= {Bi- 


VL + aDx} 


£1,20 


= {D 2 - 


-Pi4 


Pi + 


aP 5 } 


£l,38 


= {P 3 4 


-P5 4- 


eP 2 } 


A,3 = 


= {Bi- 


1-L + aPx + <5P 5 } 
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= {D 2 - 


h£>i 4 


B 5 } 




£l,39 


= {P 3 4 


-P 5 4- 


eP 5 } 


£1,4 = 


= {Bi- 


1- Di} 


£1,22 


= {D 2 - 


\- D\ 4 


P 3 + 


aP 5 } 


£l,40 


= {P 3 4 


-6P 2 } 




-£1,5 = 


= {Bi- 




£1,23 


= {D 2 - 


VDx 4 


B 5 + 


8P4} 


£l,41 


= {P 3 4 


-<5P 2 H 


-eP 4 } 


£1,6 = 


= {Bi- 


1- B 5 } 


£l,24 


= {D 2 - 


(-£>! 4 


B 5 + 


aP 4 + SP 5 } 


£l,42 


= {P 3 4 


AP 3 } 




£1,7 = 


= {Bi- 


\-B 5 + aP 4 + eP 5 } 


£l,25 


= \B 1 - 


-B 2 } 






£l,43 


= {P 3 4 


-AP 3 - 


t-aP 4 } 


£1,8 = 


= {Bi- 


H B b + eP 4 } 


£l,26 


= {B 1 - 


-B 2 + 


SL + 


aD x } 


£l,44 


-{P 3 4 


AP 4 } 




£1,9 = 


= {Bi- 


^P 3 + P 5 } 


£l,27 


= {Si- 


-B 2 + 


XL + 


aDi+8P b } 


£l,45 


= {P 3 4 


P 5 } 




-£1,10 


= {B 1 


+ P 5 } 


£l,28 


= {B 1 - 


-B 2 + 


APi} 




£l,46 


= {B 5 } 






£1,11 


= {D 2 } 


£l,29 


= {Pi- 


-B 2 + 


APi - 


f SP 5 } 


£l,47 


= {P 5 4 


-Pi + 


aP 2 } 


£1,12 


= {D 2 


+ XL + a£>i} 


£l,30 


= {B 1 - 


-B 2 + 


AP 5 } 




£l,48 


= {P 5 4 


-P 2 + 


aP 4 } 


•£1,13 


= {D 2 


+ XL + aDi + SP 5 } 


£l,31 


= {L + 


aDi} 






£l,49 


= {P 5 4 


-P5} 




£1,14 


= {D 2 


+ APi} 


£l,32 


= {L + 


aD x + 


6P5} 




£l,50 


= {Pi} 






•£1,15 


= {D 2 


+ XD 1 + SP 5 } 


£l,33 


= {Di} 








£l,51 


= {Pi4 


Ps} 




•£l,16 


= {D 2 


+ SP 5 } 


£l,34 


= {D l - 


UP 5 } 






£l,52 


= {^3} 






•£1,17 


= {D 2 


-P1+P3} 


£l,34 


= {B 3 } 








£l,53 


= W 






£l,18 


= {D 2 


- D 1 +B 3 + SP 2 } 


£l,36 


= {P 3 H 


B 5 } 















Table 3. List of 1-dimensional representative subalgebras of L. The parameters are 
e = ±1 and a, X, 5, e R, where <5 7^ 0, A > 0. 
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2.2. Classification into conjugacy classes of the subalgebra S. 

For the sake of classification, we decompose the seven-dimensional subalgebra S into 
the direct sum Q © {-Di}, with Q = {Bi, D2, B2, K, L, P5}, which is further decomposed 
as follows: 

g = A> M, 

where A is the simple algebra of the previous subsection and M. = {K,L,P 5 } is the 
ideal. Applying the method [TT1 [L2l H3] . we proceed to classify all subalgebras of £ 
into conjugacy classes under the action of the automorphisms generated by G and the 
discrete transformations fTTOj) . In practice, we can classify the subalgebras under the 
automorphisms generated by G and decrease the range of the parameters that appear 
in the representative subalgebra of a class using the Lie algebra automorphisms ( TTTj) . 
The classification results for S are shown in table H] for subalgebras of codimension 1 
and in table [5] for subalgebras of codimension 2. 



= 






Si, 14 


= |Si- 


B 2 ■+ 


-6L] 


Si, 2 7 


= W 2 - 




62D1 } 


•Si, 2 = 


{Bi + K} 




Si, 15 


= {Si- 


S 2 + 


- L + eP 5 } 


Si, 28 


= {D 2 - 


f SP 5 } 




"Si, 3 = 


{Bi + K + SDx} 




Si, 16 


= {s x - 


s 2 + 


-L + eiP 5 + ( 2 D X } 


Si, 29 


= {D 2 - 


f <5iL4 


6M 


<5l,4 = 


{Bi + L} 




Si, 17 


= {s x - 


S 2 + 


-SiL + S 2 Si} 


Si, 30 


= {D 2 - 


f SDi} 




"Si, 5 = 


{Bi + L + eP 5 } 




Si, 18 


= {Si- 


S 2 + 


-SP 5 } 


Si, 31 


= {K} 






"Si, 6 = 


{B 1 +L + ei-Di + 




Si, 19 


= {Si- 


S 2 + 


-S 1 P 5 + S 2 D 1 } 


Si, 32 


= {# + 


SD ± } 




"Si, 7 = 


{Bi +L + 5D!} 




Si, 20 


= {Si- 


B 2 ■+ 




Si, 33 


= {L} 






"1,8 = 


{Bi + P 5 } 




Si, 21 


= {D 2 } 






Si, 34 


= {L + 


SD,} 




Si, 9 = 


{Br + P 5 + 5D,} 




Si, 22 


= {^2 + 


5K} 




Si, 35 


= {L + 


eP 5 } 




Si, 10 = 


= {Si+Si} 




Si, 23 


= {D 2 + 


SiK 


+ SM 


Si, 36 


= {L + 


eiP 5 + 




<5i,n = 


= {Bi - B 2 } 




Si, 24 


= {D 2 + 


5L} 




Si, 37 


= {P 5 } 






5l,l2 = 


-- {Si -B 2 + 5K} 




Si, 25 


= {D 2 + 


L + 


eiP 5 } 


Si, 38 


= {s 5 h 


-6Dx} 




Si, 13 = 


= {Si - S 2 + 


+ <5 2 Si} 


Si, 26 


= {D 2 + 


L + 


eiP 5 + 


Si, 39 


= {Di} 







Table 4. List of 1-dimensional representative subalgebras of S. The parameters are 
e, ei, e 2 = ±1 and 6, 5\ 1 5 2 € R, <5, #i,£ 2 7^ 0. 



3. Invariant and partially invariant solutions. 

Since equations ((Ua) and (JUb) do not involve the velocity components u and t> , they 
can first be solved for 9 and a. Next, the result is introduced into the system formed 
by equations ([Be) and ([Ud) and we look for the solution of this system for the velocity 
components u and v. This system always admits the particular solution 

u = hy + b 2 , v = —b 1 x + b 3 , 6i,6 2 ,6 3 eR, (15) 

obtained by requiring that the coefficients of the trigonometric functions in ([He) vanish. 
The velocity field defined by (|15|) forms concentric circles. Therefore, this solution does 
not establish any relation between the flow velocity and the strain involved in the plastic 
material. Therefore, it is not an interesting result by itself from the physical point of 
view. Nevertheless, since the PDEs (JTJc) and ([Ud) are linear (assuming that 9 is known), 
they admit a linear superposition principle. Then, we can add the solution (I15j) to any 
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^2,32 


= Pi - Bo 
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+ dK,Di\ 
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^2,2 " 
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^2,33 


= -TBi - Bo 

— J --'2 
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— i Pi - Bo 


+ L + 




*^2,4 _ 


- i Pi Do 4- e 1 L 4- f 9 P* 1 


^2,35 


— j R 1 _ R 9 


+ <SP 5 ,Bi} 
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= iBi Bo + 

\ J -- / li J - y 2 1 u ± of 


^2,36 


= - Bo 


+ SDi,K + aD 1 } 


<^2,6 ~ 


= iBi D? + SiK + <5oBi i 


*^2,37 


= {Bi - B 2 


+ SD 1 ,L + aiD 1 } 


^2,7 ~ 


= iPi Do 4- eiL 4- eoP^l 


^2,38 


= {Bi - B 2 


+ SD 1 ,L + e 1 P 5 + e 2 Bi} 


c 

t - , 2,8 " 


- 1 Pi Do -4- fi L 4- fo Pc 4- Di \ 


c 

^2,39 


= {B 1 - B 2 


+ SD 1 


,P5 + aiB!} 


c 

t - 7 2,9 " 


= iPi Do 4- SiL 4- <5oDii 

\ J --'l 1 ^2 1 U l- L; t u 2 ^ 1 j 


c 

t - , 2,40 


= {D 2 ,K} 






^2,10 


= -f Pi Do 4 f 1 4- FO Pr -1- Dl 1 
\ J --'l; -^2 I c l iJ I C 2 J D ^ C 3- Ly l J 


c 

( - , 2,41 


= {D 2 ,L} 






^2,11 


\ J --'li- Ly 2 I u i- 1 > u 2 J - y L f 


^2,42 


= {D 2 ,L + 






^2,12 


= \B^ Bo + 5Di\ 


^2,43 


= {D 2 ,P 5 } 






c 

t - , 2,13 


— 1 1 ; J 


^2,44 


= \d 2 ,da 






^2,14 


— / Ri r\ 

— \ 1 ? -^x 


^2,45 


= {Bs + ^Bx} 




^2,15 


= IBi L + eP*\ 


*^2,46 


= {D 2 + 5L,D 1 } 




^2,16 




^2,47 


= {Ba+i + ePs,!?!} 


^2,17 


— / r, n, i. 


^2,48 


= {B 2 + «5P 5 ,B 1 } 




^2,18 




<^2,49 


= {B 2 + ABi,if + aBi} 


C 

t - 7 2,19 


= iBi + L Bii 


C 

^2,50 


= {B 2 + AB 1 ,i + aB 1 } 


^2,20 


= -TBi + L + eA Do + 2eoK\ 


^2,51 


= {B 2 + AB 1 ,i + eiP 5 + e 2 B 1 } 


c 

t - , 2,21 


— f r 1 4. t 4. frP, nj 

— \ ^ 1 T T ci 5, .L^l f 


c 

t - , 2,52 


= {B 2 + AB!,P 5 - 


haBi} 


^2,22 


= iPl 4- P^ Dl 1 


^2,53 


= {L,P h + aD 1 } 




(S2 23 


= {Bi + Di,K + aDi} 


(S2 54 


= {L + 6D 1 ,P 5 } 




^2,24 


= {B 1 +D 1 ,L + aD l } 


^2,55 


= {L + e l D 1 ,P 5 + 




^2,25 


= {B 1 +D u L + e 1 P 5 + e 2 D 1 } 


^2,56 


= {K,L + eP 5 } 




<^2, 26 


= {B l +D 1 ,Ps + aD 1 } 


^2,57 


= {L,P 5 } 






^2,27 


= {B 1 -B 2 ,K} 


£2, 58 


= {^,^1} 






^2,28 


= {B X -B^L) 


£2,59 


= {L,D 1 } 






£2,29 


= { J B 1 -B 2 ,L + eP 5 } 


<^2, 60 


= {L + eP 5 ,D 1 } 




^2,30 


= {B 1 -B 2 ,P 5 } 


<$2,61 


= {P 5 ,D 1 } 






£2,31 


= {B 1 -B 2 ,D 1 } 











Table 5. List of 2-dimensional representative subalgebras of S. The parameters are 
e, ei, e 2 = ±1 and a, A, 5, 5\,5 2 e M, 6, Si, S 2 7^ 0, A > 0. 

solution, for the velocity components, of the system ([He), ([Ud), corresponding to given 
solutions 9 and a of the system ([Ua), ([Ub). Consequently, we can satisfy a much broader 
family of boundary conditions. 

3.1. Symmetry reduction for the representative subalgebra B\ 

Consider for illustration the one- dimensional representative subalgebra generated by the 
infinitesimal generator 

Bi = —vd x + ud y . (16) 

Since no derivative with respect to variables u,v,9,a appears in B\, it follows that 
these variables are all invariants of the subalgebra generated by B\. In order to obtain a 
complete set of functionally independent invariants, one can include also the symmetry 
variable 

£ = ux — vy. (17) 

We look for a solution of the form 



(18) 
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where £ is defined by (ITT]) . Replacing f llSp into the original system f )16p and assuming 
that 1 - xF'(g) - yG'iX) ^ 0, where F'(£) = dF(£)/d£, etc., so that we can use the 
Implicit Function Theorem, we obtain a system of equations with a reduced number of 
independent variables, where the functions F, G, T, S are to be determined. This system 
takes the form of four coupled ODEs: 

F(£)S'(£) - (cos(2T(£))F(£) + sin(2T(£))G(£)) T'(£) = 0, 

G(£)S'(£) - (sin(2T(£))F(£) - cos(2T(£))G(£)) HO = 0, 

(G(£)F'(£) + F(0G'(0) sin(2T(£)) + (F(£)F'(£) - G(£)G"(£)) cos(2T(£)) = 0, 

F(£)F'(£) + G(£)G'(£) = 0. 

The solution of this system is 

F(£) = Cl cos(T(£)), G(£)= Cl sin(T(£)), S(£) = T(£) + c 2 , (20) 

where T(£) is an arbitrary function of a single variable and C\,C2 are constants of 
integration. The solution is obtained by replacing expressions (120]) for F, G, T, S into 
equation ( TT8]) . The solution of the original system is given implicitly by relations 

u = Ci cos(wx + yv ), f = Ci sin(wa; + yt> ), (21) 

while the angle 9 and the pressure a are defined by the choice of the arbitrary function 
T as follows: 

9 = T(ux + yv), a = T(ux + yv) + C2- (22) 

Since by defining 6 through a certain choice of T, we also determine a, it follows that 
equation ( 1221) is a relation defining the pressure a in terms of the angle 9 or vice- 
versa. Moreover, we can see from equation (12T|) that the sum of squares of the velocity 
components u and v, is constant. Therefore, since the material is incompressible, the 
solution preserves the kinetic energy of the plastic material, i.e. 

2,2 2 
U + V = c x . 

For the purpose of illustration, consider the function: 

T(0 = (l/2)arcsin(£), (23) 

This particular choice of T allows us to solve relations (12TT) in order to find the velocities 
u and v explicitly in terms of x and y. The obtained formulas can be expressed in terms 
of radicals and are very involved. Therefore, they are omitted here. Nevertheless, these 
formulas can be used to trace the vector fields corresponding to solution ( I2T|) . where T 
is defined by f l23l) . An example of such a tool is presented in figured] for a flow velocity 
Ci = 5. The feeding velocity used is (Uq, Vq) = (4.30,2.55) and the extraction speed is 
(Ui, V\) = (—4.30, 2.55). The boundaries of the extrusion die are chosen in such a way 
that they coincide with the flow lines of the velocity field. Therefore, they are solutions 
of the equation 

dy/dx = v(x,y)/u(x,y). 
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For figure 1, the inner boundary corresponds to the initial value (xq, yo) = (—0.5, —0.35) 
and the outer boundary to the initial value (xo,yo) = (—0.43,-0.46). The curves C\ 
and C2 are the limits of the plasticity region with respect to the entrance and exit of the 
extrusion die. They are solutions of equation (J2j), where, (Uq, Vq) is replaced by (U\, V\) 
for C2. In order to define the limit of the plasticity region at the ends of the boundary of 
the tool, the initial data used to trace the curve C\ are (x , yo) = (—0.5, —0.35) while for 
C2 they are (xo,yo) = (—0.5,0.35). Numerical integration has been used to identify the 
boundary of the tool and the limits of the region of plasticity. This type of extrusion die 
can be used to bend a rectangular rod or a slab of a ideal plastic material. The average 
pressure and the angle 6, which define the strain tensor inside the tool, are evaluated 
by formulas (122]) . where T is defined by (l23l) . 




Figure 1. Extrusion die corresponding to the solution (f^Tj) . 



3.2. Symmetry reduction for the representative subalgebra K 

As an example, we find a partially invariant solution corresponding to the subalgebra 
generated by generator K which admits the following invariants: 

f = xycos(2fl) - (l/2)(x 2 -?/ 2 )sin(2fl), F = xu + yv. (24) 

and 

S = 2 - a 2 , G = uvcos(26) - (1/2) (u 2 - v 2 ) sin(20). 

In order to obtain a PIS, we use only the two invariants given by (124)) . We begin by 
inverting the first relation in (124"1) in order to find 9 as a function of £. Next, we introduce 
the result in the first two equations of system ([T]). Then, comparing the values of the 
mixed derivatives of a(x, y) with respect to x and y, we obtain the following PDE for 
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the quantity £ 

(£« - U (((* 2 + y 2 ) 2 - ^em* 2 - y 2 ) - xvVW+W^)) 

Axyj + (x 2 - y 2 ) y/(x 2 + y 2 ) 2 - 4£ 2 2 2 2 

xyy [x + y ) — 4£ — [x z — y z )^ 
x ((x - y)£ x + (x + y)f„) - 4(x 2 + y 2 ) 2 ^ (x£ x +y£ v -£) = Q 

There are two particular solutions for £ to equation ( |25|) defined by 

£ = i(x 2 + ey 2 ), e = ±l. (26) 



Let us consider the case e = 1 and introduce this solution into the first relation ( 1241) . 
Solving for 9, we obtain: 

1 / ' x 2 — y 2 * 



= — arctan — . (27) 

2 \ xy ) K > 

The mean pressure a is found by quadrature from the first two equations (JTJ) in which 
we have introduced the solution (1271) for 9. The result for a is: 

a = -(l/2)\n(x 2 + y 2 )+c u (28) 

where c\ is a real integration constant. 

Using the form of the second invariant in (|24|) . we look for a solution for the 
components u and v of the velocity, with the form 

u = (y/x)v-F(0, (29) 

where the symmetry variable £ is given by (126]) . By replacing 9 given by (j27|) and 

u by (129]) into the system composed of the last two equations in ([T]), then using the 

compatibility condition of the mixed derivatives of v with respect to x and y, we obtain 

the condition that F = C2, where C2 is a real constant. The solution for u and v is then: 

C2X Ciy 
u = — J + c 3 y + c 4 , u = — =■ - C 3 X + c 5 , (30) 
x + y x + 2/ 

where C3 . . . , C5 are real constants of integration. Note that in the case when 04 = 05 = 
and C2 7^ 7^ C3, the flow lines form logarithmic spirals centered at the origin. 

An example of velocity fields is given in Figure [2] for parameters C2 = —1, C3 = —2, 
C4 = 4 and C5 = 1 for solution ( 130]) . The chosen region, [—1, 1] x [—1, 1], includes the 
singularity at the origin. Corresponding to this solution for the same parameters, an 
extrusion tool is given in Figure [3] for the feeding and extraction velocities (Uq,Vo) = 
(5.5,0) and (Z7l,Vi) = (3,3) respectively. The curve C\ is the limit of the plasticity 
region at the entrance of the extrusion die and C2 has the same significance at the exit 
of the extrusion die. The upper contour of the extrusion die is a solution of dy/dx = v/u, 
with u and v defined by (130]) . for an initial value y{— 0.5) = —0.8 while, for the lower 
contour, we have used the initial value y{— 0.7) = —0.95. 



Group analysis of an ideal plasticity model 



13 



Sift 

■% "Si 



"** 



*** ^ 
^ y y 
»* y 



-~4 



^ \ \ 



■«* 



/ - 



/ / 
/ / 



/ / 

/ / 

/ / 

/ / 

s / 

/ / 

/ / 

/ / 

/ / 

/ / 



/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 
/ 

/ X 



■ y 

■ s" 

■ y 

■ y 



\ \ 

t t 

/ / 

S / 

s s 

y s 

y y 

y y 

y y 



I / A 

/ / / 

/ / / 

/ / / 

/• s s 

y s s 

y y y 

y y y 

y y y 



/ / 

/ / 

/ / 

/ / 

/• S 

y y 

y y 

y y 

y y 



/ / 

/ / 

/ / 

s / 

s / 

s s 

y y 

y y 

y y 



Figure 2. Extrusion die corresponding to the solution (f5U|) . 




Figure 3. Extrusion die corresponding to the solution (|30|) . 



3. 3. Similarity solution for the angle 9 and corresponding pressure a 

In this section, we find solutions of the system ([1]) for which the angle 9 is a similarity 
solution. We propose the solution for 9 in the form 

e(x,y) = J(£(x,y)) i (31) 

where the symmetry variable is of the form 

Z{x,y)=y/x. (32) 

One should note that the solutions, obtained by assuming the hypotheses (}3~Tj) and (|32|) 
on their form, are more general than the invariant solutions corresponding to subalgebras 
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which admit y/x as a symmetry variable. This is so because the invariance requirement, 
for a given subalgebra, leads to constraints on the form of a, which is not the case here. 
The introduction of ( 13T|) . with £ defined by ( 1321) . in the system (JTJa), flUb), leads to the 
system 

a x (x, y) = 2(k/x) (-£(z, y) cos( J (f (x, y))) + sin(£(x, y))) J'(£(x, y)), 
^(z, y) = -2(fc/x) y) sin( J(£(x, y))) + cos( J(£(z, y)))) J'(£(z, y)). 

Considering the compatibility condition on mixed derivatives of o relative to x and y, 
we deduce from (1331) the following ODE for the function J: 

((e - 1) sin(2J(0) + 2£cos(2J(0)) J"(0 + 2 (£sin(2J(0) + cos(J(0)) -AO ^ 
+ 2 (-2 sin(2 J(0) + (£ 2 - 1) cos(2 J(0)) ( J'(0) 2 = 0, 

which has the first integral 

((e - 1) sin(2 J(0) + 2£ cos(2 J(0)) J'(0 = c 1? (35) 
where Ci is an integration constant. There are two cases to consider to solve the equation 

(El. 



zj. If Ci 7^ 0, the solution of ( I35j) is given in implicit form by 



" « - tan V^ 1 '^ - J «» I = 0, (36) 



(tan(J(e))e + l)ft-{ + tan(J({)) ^ Cl 

where C2 is an integration constant. The solution for a is obtained by integrating the 
system (133]) and taking into account the first integral (135]) . We find 

a(x, y) = k {£{x, y) cos(2 J(£(x, y))) - sin(2 J(£(x, y))) - 2 Cl ln(x)) + c 3 , c 3 eR (37) 

where y) is defined by ( 132]) and J(£) is a solution of ( 136]) . So, 9(x,y) = J(y/x) and 
cr(x, y) are given by (137]) and are solutions of the system (JTJa), ([Ub), if the function J(£) 
satisfies the algebraic equation ( 136]) . 



ii). If the constant ci = in equation (13"6"]) . then the solution of (131]) for J is 



J(0 = -^arctanf_^L). ( 38 ) 



We subsequently solve (133]) for a{x,y) considering (138]) . which leads to the solution of 
the system fllJa), jUb), given by 

9(x, y) — - arctan I — I , a(x, y) = —2k arctan ( — ) + ci- (39) 

2 \x l — y l ) \xJ 
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3.3.1. Additive separation for the velocities when c\ 7^ 0. Already knowing the solution 
9(x, y), a(x, y), in the case where C\ 7^ 0, we still have to compute the solution for u and 
v. One way to proceed is to suppose that the solution is in the additive separated form 

u(x,y) = f(x,y) + F(£(x,y)), v(x, y) = g(x, y) + G(£(x, y)), (40) 

where £(x,y) = y/x. We introduce ([4*0]) into the system ((TJc), dUd), which gives 

(sin(2 J(0) - cos(2 J(0)) tnt) - (cos(2 J(0) + sin(2 J(£))) £G"(0 
+ (Uy + 9x) sin(2 J(0) + {f x - g y ) cob(2 J(0)) z = 0, 

(/x + ^)x + G'(O-^'(O=0, (42) 

We must determine which functions / and g will reduce equations (HTj) . (|4"2|) . to a system 
of ODEs for the single- variable functions F(£) and G(£). To reach this goal, we first 
use as annihilator the infinitesimal generator (1/2) (Di + D 2 ) defined by (T4j) that we 
apply to the equations (HI]) . (l4*2|) . to eliminate the presence of the functions F and G. 
Indeed, the operator (l/2)(Z?i + D 2 ) annihilates any function of £ = y/x. So, we obtain 
as differential consequences some conditions on the functions / and g. We can assume 
that f x (x,y) + g y (x,y) 7^ 0, otherwise we can show that the only possible solution is 
the trivial constant solution for u and v. Under this hypothesis, the previous conditions 
read 

f x = -9y + Ui)x-\ (43) 

fy = -9 X +9yU0 + U0x-\ (44) 

(g y + (A + D 2 ) (g y )) (C 2 (0 sin(2 J(0) - 2 cos(2 J(0)) = 0, (45) 

where the functions of one variable Q, i = 1, 2, 3 are arbitrary. Since the left member of 
151 is composed of two factors, we must consider two possibilities. 



(a). We first suppose that 

g y + (D, + D 2 ) (g y ) = 0. (46) 
In this case, we find that the functions / and g take the form 



f{x,y)=-J 7 d£ + u 4 \n(y) - u x y + w 5 , 



rt(*,y) Cl (0-C(0 

(47) 

g(x, V) =C4(0 + ^2 ln(a;) + w x x + w 3 , 



where the functions Ci(0> Ci(0 are arbitrary and the functions C2 (C) 5 Cs(0> were chosen 
to solve the compatibility conditions on mixed derivatives of / relative to x and y. We 
now introduce the solution ( 1471 in the system ( l43l . ( 1441) . which leads to an ODE system 
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for F and G, that we omit due to its complexity, and for which the solutions take the 
form of quadratures 



F(f)- I f C'(D I 0(0 I (^4 + ^0sin(2J(Q) 

M ^ )+ £ + m 2 -l)sin(2J(0) + 2£cos(2J(0)) ^ + C5 ' 



J (-Ci(0 + ^'(0)^ + c 4 , c 4 ,c 5 G 



(48) 



The last step is to introduce (|47p and (|48|) in the Ansatz (I4T)]) . Then, the velocities are 

l c 5 cos(2J(y/x)) /•»/« c 6 sm(2J(0)J'(0 ^ , 
2 ci J ci4 

cr 1 c 5 esin(2J(0)J'(0^- (2c 1 )- 1 (c 6 cos(2J(0)) + c 8 , 

(49) 

where the c, are integration constants. So, we obtain a solution of the system ([T|) by 
defining the angle 9 by (j31jl . (1361) . the mean pressure a by ( 1371) and the velocities w and 
v by (jUJ), with £(x,y) = y/x. 

(b). Suppose now that the condition (I4"5j) is satisfied by requiring 

C 2 (0 = 2cot(2J(0). (50) 

Then, applying the compatibility condition on mixed derivative of / relative to x and 
y to the equations (143]) . (144|) . and considering £ 2 given by (150|) . we conclude that the 
function g must solve the equation 

g xx (x, y) + 2 cot(2 J(0 W*. 2/) ~ 0w(*> 2/) + 4s_1 ^ (f + cot 2 (2 J(0)) J'(Z)g v (x, y) 

+ x- 2 (cm + cm + cs(0) = o- 

(51) 

It's a hyperbolic equation everywhere in the domain where J is defined. So, we introduce 
the change of variable 



sin(2J(0) 

(X ' y) =X 6XP 1 / l + cos(2J(0)+^in(2J(0) ^ 

, sin(2J(Q) 

=xexp / — — 

-l + cos(2J(f)) + fsin(2J(0) 



(52) 



which brings the equation ( 1511) to the simplified form 
ci / sin(2J(0,^))^ sin(2J(0,^))^ 



2 V^(cos(2 J(0, V))) + 1 0(cos(2 J(0, VO) - 1) 

1 sin(2J(</>, ^))£ 2 - sin(2 J(0, + 2£ cos(2 J(0, jj) (C 3 (Q + C( (0 + ^(0) 

4 J(0,V) 

(53) 



0. 
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where J(0, %p) is defined by 

J =(ci/4)(<A -</>)• (54) 
To solve the equation (1531) more easily, we define the function £3 by 

c 3 (0 = r 1 (-Ci(0+^(0+wi). (55) 

So, the solution of (|53|) is 

s(0,V) = -(1/2) (w a - (1/2) ln(^/0)) cos ( Cl ln(^)) -(1/4)0^ sin ( Cl ln(^/0)), (56) 

which, returning to the initial variables, takes the form 

g(x,y) = -(1/2) ( Wl - (l/2)cr 1 J(y/x)) cos(2 J(y/x)) - (1/4) of 1 sin(2J(y/x)) + w 2 . 

(57) 

After the introduction of the solution (1571) for g, the function / is given by quadrature 
from the equations (|43p . (jUJ). The obtained solution for / is 



f(x,y)=\ I — )~J — ^— dC+(ci+l)wiln(2/)+W3- 

(58) 

We now introduce (157j) . (1581) . in (14T1) . (14*21) . and get F and G by quadrature in the form 

F(n f ax ' y) Ci(£) , (^ + J(0)sin(2J(0)^(0 

F(0 = i + i ^ (59) 

G(0 = -(l^cf 1 ((wj + J(0) cos(2J(0) + 2 sin(2J(0)) . 



Finally, the substitution of (l57j) . (|58|) and (159]) in (14D|) provides the solution to ([He), 
(ffld): 

/ / ^ WN r /a; C4(ci + l)sin(2J(0)J'(0^ 
«(x, y) =(ci + l)c 4 ln(y) + / ^ \ ^ dg + c 5 , 

J Cl4 (60) 

c 4 (ci + 1) cos(2J(y/x)) 
*>(z, y) = + c 6 , 

Cl 

where the Cj are integration constants. So, we have a solution of the system ([1]) by 
implicitly defining the angle 9 by ( 13"T1) . ( 1361) . the mean pressure by cr by ( 1371) and the 
velocities w and v by (160|) . with £,(x,y) = y/x. 

3.3.2. Additive separation for the velocities u and v when c\ = 0. Now, we consider 
the case where c\ = in (j35p . Then the solutions for 9 and cr are given by (13"9"j) . We still 
suppose that the solution for u and t> is in the form (1401) . The procedure is the same as 
for the previous case until we obtain the conditions ( l43l) . (l44j) and ( 14"5l) . We must again 
consider two distinct cases. 
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(a.) We first suppose that the condition (146 D is satisfied. Then, the functions / and g 
are defined by 

f(x, y) = f iX ' y \uO + - Ca(£R - my + w&, (gl) 

9(x, y) =(4,(£(x, y)) + ujix + u 2 , i = 1, 2, 3, 

where the £j, = 1,2,3, are arbitrary functions of one variable and, to simplify the 
expression for / and g, we have chosen Ci(0 = (1 - (2(^)6 ~ £ 2 ) C*(0 ~ fCs(0- We 
substitute (l6"Tj) in equations (l4"Tj) and (1421) to determine F and G. We conclude that 
F(£) is an arbitrary function, while G(£) is expressed as the quadrature 

G(0 = J (F(Oe + (foCfl - 1 + e 2 ) CI - eCs(O) df. (62) 

We finally obtain the solution m and u by introducing ( 161~|) . ( 1621) in (HDIl and, to simplify, 
by choosing 

C 3 (0 = -(C 2 (0 + OC(0, 

which gives 

u(x,y) = - c 2 y + F'(y/x), 

v(x,y) =c 2 x + (y/x)F'(y/x) - F(y/x), 

where F is an arbitrary function of one variable. A solution of the system ([1]) consists 
of the angle 6 and the pressure a defined by ( 1391) with the velocities defined by ( 1631) . 
For example, if we choose the arbitrary function to be an elliptic function, that is 

F(£) = cn ((1 + cosh(arctan(6 2 0)) _1 ,q), < p 2 < 1, 

and we set the parameters as h\ = Air, c 2 — 0, p — 1/2, then we can trace (see figure 
H|) an extrusion die for a feeding speed of (Uq, Vq) = (0, —0.94) and an extraction speed 
(Ui,Vi) = (0,-0.94). The curve C\ on the figure H] delimits the plasticity region at 
the mouth of the tool, while the rc-axis does the same for the output of the tool. This 
type of tool could be used to undulate a plate. We can shape the tool by varying the 
parameters. For example, we can spread the bump by decreasing the parameter h\. 
Moreover, one should note that if the modulus g of the elliptic function is such that 
< g 2 < 1, then the solution has one purely real and one purely imaginary period. For 
a real argument Xi we have the relations 

-1 < cn(x, g) < 1. 



(b). Suppose that Ca(0 is defined by (150]) and for simplification we choose in particular 

Cs(0 = (64) 
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Figure 4. Extrusion die corresponding to the solution (|39|) . (|63|) . 



Applying the mixed derivatives compatibility condition of / to the equations 
we get the following ODE for the function g: 

g xx -g yy -2 cot(2J(0)g xy - 4 ^f (0 ^ = 0. (65) 

ism (2J(£)J 

By the change of variable 

€(x,y)=y/x, v(x, y) = x 2 + y 2 , (66) 



we reduce the PDE ( l65i) in term of x and y, to the much simpler PDE in term of £ and 
V, 

9*, + 4^7 = °> (67) 



which has the solution 

g^v) = U0 + ^==i (68) 

where £4 and £5 are arbitrary functions of one variable. Then, we find the solution for 
/ by integration of the PDE (j4"5|) . (jUJ), with £ 3 given by (J6i 



f{x,y) = - 7 , . = +c 2 . (69) 

By the substitution of (1681) and (I6"9"j) in the equations (141 p . (j42p . we find that F is an 
arbitrary function of one variable and G is defined by 

G(t{x,y)) = f iX,V) (-Ci(0 + F'(00^ + c 3 . (70) 

So, we introduce (|68|) . (I6"9"j) and (I7D|) in (|4U|) and after an appropriate redefinition of £1, 
(4 and (5, the solution of (EDc), flUd) is provided by 

= K'(y/x)-yH(x 2 +y 2 ) + c 2 , v(x,y) = -K'(y/x)+£K(£)+xH(x 2 +y 2 ) + c 1 , 

(71) 

where H, K are arbitrary functions of one variable. The velocities ( 17TT) . together with 
the angle and pressure defined by (1391) . solve the system ([1]). For example, a tool 
corresponding to the solution ( 17TT) with #(77) = 2 exp(— O.I77) , K(£) = £ and for feeding 
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and extraction speed given respectively by (Uq,Vq) = (1.05,0) and (Ui,Vi) = (1.05,0). 
It is shown in figure [5j The plasticity region limits correspond to the curves C\ and Cg. 
This tool is symmetric under the reflection x \- > —x. Moreover, the top contour of the 
tool almost makes a complete loop, and this lets one suppose that we could make a ring 
in a material by extrusion. 
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Figure 5. Extrusion die corresponding to the solution (|71[>. 



3.3.3. Multiplicative separation for the velocities u and v when c\ ^ 0. Consider the 
solutions of (jUa), (JT]b) , given by the angle 8(x,y) = J(y/x) with J defined by (l36l) 
and the pressure a defined by (I37I) . We require that the solutions for u and v be in the 
multiplicative separated form 

u(x, y) = f(x, y)F(£(x, y)), v(x, y) = g(x, y)G(£(x } y)), (72) 

where £(x,y) = y/x and /, g, F, G are to be determined. The Ansatz (1721) on the 
velocities brings the system ([He), ([Ud) to the form 

[f v F + x-'fF'iO + g x G(0 ~ x-^gG'®] sin(2 J(£)) 

+ [f x F(£) - x- x UF' - g y G - x-'gG'] cos(2J(0) = 0, (73) 

[f x F($ + g y G(0) - CfF'(0 + gG'(0 = 0. 



To reduce the PDE system f[73|) to an ODE system involving F and G in terms of £, we 
act with the operator (l/2)(Z?i + -D2), defined by ([!]), and annihilate the function of £ 
present in (173|) . This leads to conditions on / and g that do not involve F, G and their 
derivatives. There are three cases to consider, that is 

(o , (1/2X* + ft) {9M + 0, (ft + ft) (fig^ggg)) * 0, 

(6) (l/2)(D 1 + D 2 )(x(/ 1 ,/^) = 0, (74) 

(0 (1 / 2 )(ft + ft)„v s) ,o, 

In this paper, we present the details for the cases (a) and (b). 
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(a). We first suppose that the condition (1731 a) is satisfied. In this case, the function 
/ must be a solution of the PDE system 

/. = - s-^CiCfl/ + - x ~ l (UOCiiO - C 2 (0) 9, 

f y =z -1 Ci(0/ + C 2 (0fc - 2 C 2 (0 cot(2j(0K (75) 

+ (C 2 (OCx(0 - C 2 (0) xsin(2J(0) 5, 



where £ 2 are two arbitrary functions of one variable. For the system (175]) to be 
compatible, the function g must satisfy the PDE 

g xx - 2 cot(2 J(£))g xy - g yy + 2 X - 1 C 2 (0~ 1 (C + cot (2 J(0)) Cs(0^ 
+ 2X- 1 [(£ - cot(2J(0)) Ca(0 - 2^sin- 2 (J(0) 0„ 

+ x- 2 C 2 (0- 1 ([4^'(0sin- 2 (2J(0) + (e 2 - 1 + 2^cot(2J(0)) C'i(0 (?6) 
- 2£ - 2 cot(2J(0)]C 3 (0 - (e - 1 + 2Ccot(2 J(0)) Cs(0)ff = 0, 



where we used the notation Cs(0 = C2 (C) Ci (0 — C 2 (0 to shorten the expression. The 
equation (176]) is difficult to solve for arbitrary (1, ( 2 , but if we make the particular choice 

C 2 (£)=u; 2 exp(Ci(0), u 2 ER, (77) 

then the PDE ( 176]) reduces to 

g xx - 2 cot(2 J(Z))g xy - g yy - 4x~^ J'(0 sin" 2 ( J(0)^ y = 0, (78) 

which is solved by the function 

<7(x, y) = U3X + U4 cos(2J(y/x)) + W5, W3, W4, G K. (79) 

With g given by (1791) and C 2 (0 by (ITTj) . the system (175]) is compatible. Consequently, / 
is expressed in term of a quadrature. We find 

f(x,y) =2exp(( 1 (y/x))L0 2 L0 4 f 1 sin(2J(£)) J'(£R (gQ) 

+ exp(Ci(y/x)) (w 2 w 3 + 2w 2 W4Ci ln(y) + w 6 ) 
With / given by (IHUj) and (7 by (1791 . the solutions for F and G of the system (173|) are 

F(0=w 1 exp(-Ci(0), G(0=w 2 . (81) 

By introducing (179 p . (18"U|) . (18 ip in (172]) and redefining the free parameters cjj, i = 1,2, 3, 4, 
the solution of flTJc), dUd) for the velocities u and t> is 

/j//as 
r 1 sin(2J(0) /(O^ + c 5 2/ + 2c lC4 In(y) + c 6 , 

f (x, y) = — C5X — C4 cos(2J(y/x)) + C7, 

where the q, i = 1, . . . , 7, are integration constants and J is defined by (13"6]) . So, we 
have a solution of the system ([T|) composed of the angle 9 in the form (j3~Tl) with J given 
implicitly by (155]) together with the pressure o (J57J and the velocities 
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(b). Suppose now that the condition ()74l c) is satisfied. In this case, the solution g takes 
the form 

g(x, y) = hi(x)(i(£(x, y)) (83) 
and the function / must be a solution of 

f y = - 2x- 1 h 1 {x)Cs(0 cot(2J(0) + e-<>® ^ J"'* e*>® - u 2 ^j h[(x) (84) 

We omit, due to its complexity, the expression of the compatibility condition on the 
mixed derivative of / relative to x and y. Nevertheless, making the specific choice 

«Q - h ( _ ^l% c0B(2m ) ~ HUQ), (85) 



the system ( (841) turns out to be compatible and the solution of ( 184j) is 

f(x, y) =Uv/x) (c- 1 (1 - (y/xf - 2{y/x) cot(2J(y/x))) ["'* sin(2 J(£)) J'(£)d£ 



+ [1 - (y/x) 2 - 2(y/x) cot(2J(y/x))] (bay - u 4 ) 

(86) 

The substitution of ( 1861) , ( 1831) and ( 185|) in the system (173|) results in an ODE system for 
F and G, omitted due to its complexity, which has the solution 

F(n = msin(2J(fl) 

W Ca(0 «£ 2 - 1) sin(2J(0) + 2ecoe(J(0)) ' (87) 
G'(0=Ci(0" 1 (^icos(2J(0)-u; 2 ). 

We finally obtain a solution for the system ([He), ([Ud) by the substitution of (186]) and 
flH?]) . with fc x (ar), (2(6 > defined by (|55]l. in (1721. This leads to 

y) = -c 4 ln(y) - ^ f" 1 sin(2 J(£)) J'(£R + c 5 , 

= (l/2)c^Vcos(2J(?//a:)) + c 6 , ci,c 4 ,c 5 ,c 6 e E. 

So, the system (CQ) is solved by the angle 9 in the form ( 13T1) with J implicitly defined by 
(136|) . together with the pressure a (157)1 and the velocities (IHHj) . 



3.3.4- Multiplicative separation for the velocities u and v when C\ = 0. Consider now 
the case where C\ = in (|35|) so the solution of (JTJa), dUb), for and o is (|39|) . We 
suppose that the velocities u and v are in the form (172]) . Introducing this form for 
the velocities and # defined by ( l39i) in the equations ((Uc), ((Ud), leads to the system 
(173|) . which reduces to an ODE system for F and G if the functions / and g satisfy the 
condition (1741 . The three different constraints ( 1741) must be considered separately. 
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(a). In the first case, where we consider that the conditions (1731 a) are satisfied, the 
functions / and g must satisfy the system (1751) . (176]) . Changing the variables (x,y) to 
the new variables (£, rf) defined by 

t(v,y)=y/x, r)(x, y) = x 2 + y 2 , (89) 

and considering 9 given by ( 1391) . the system ( 1751) . ( 1761) . becomes 

h = -Ci(0/ - r'cm - r 1 (£(0 - Ci(oc 2 (0) <?, 



(90) 



— + ^® + Ci(0 + ^ 2 + ir 1 = o, (91) 

#r, C 2 (£) 

where Cl(0> C 2 (6) are arbitrary functions of one variable. The solution of the system 
(190]) . (}9Tj) . for / and (7 as functions of £ and 77 is 

J c\e + 1) [c 2 (o + (ci(o + e(e 2 + i)- x c 2 (o) ^ (o 

v (92) 

+ C 2 (0tf'(0R-"iJ, 

where the functions and -ff (77) are arbitrary functions of one variable. Introducing 

/ and g expressed in the initial variables x, y, by the substitution of (j89p in (1831) . the 
system (173]) is reduced to an ODE system for the functions F and G in term of £ which 
has the solution 

= utef^W, G(0 = -^iC 2 (0e /Cl(0d? . (93) 

Finally, redefining 

C 2 (0 = -(^(or'exp Ci(O^) Q(0, = ^VvP(v) 

where Q(£), P{f]) are arbitrary functions, and doing the substitution of (192}) and (193]) 
in (172|) . we obtain the solution of (JTJc), (JT]d) , for velocities u and v given by 

y) = yP(x 2 + y 2 ) - x^yQiyjx) Jf'* r 1 ((£ 2 + W(0 + £<?(£)) ^ + c 2 , 
v(x, y) = Q{y/x) - xP(x 2 + y 2 ), 



(94) 



This solution for velocities together with 9 and a defined by (139]) solves the initial system 

(HD. 
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(b) . Consider now that the functions / and g satisfy the constraint (1741 b). Then they 
take the form 

f(x,y) = &(y/x), g(x,y) = ( 2 (y/x), (95) 

where (j and £2 are arbitrary functions of one variable. Since the functions / and g 
depend only on the symmetry variable and the velocities u and v have the form (|72|) . it 
is equivalent to consider that 

u{x, y) = F{C(x, y)), v{x, y) = G{£(x, y)). (96) 

If we suppose that u and v are in the form ( 1961) . then the solution of (CD) consists of the 
angle 9 and the pressure a given by (139]) together with the velocities 

u(x, y) = F(y/x), v(x, y) = J" £P'(£K + c 3 , (97) 
where F is an arbitrary function of one variable. 

(c) . The third case to consider is when / and g obey the conditions (1741 c). so they 
take the form 

f(x,y) = u 3 x 1+ ^(y/xY 1+ ^ 2 (x- 1 (x 2 + y 2 )r/\ g(x, y) = x 1+ "Xi(y/x), (98) 

where Ci(0 is an arbitrary function of one variable. Then we introduce (l98l) in ( 173]) and 
solve for F and G. The solution is 

Finally, substitution of /, g, F, G, given by (198]) and (199|) . in ([75]) gives, after redefining 
the parameters in a convenient way, the solution for u and u of the equations ([He), 
dUd), 

y) = c z y(x 2 + y 2 )^/ 2 , v(x, y) = -c 4 x{x 2 + y 2 )^ 2 , (100) 

where C3, C4 are integration constants. The velocities u and v together with the angle 9 
and the pressure o given by ( 139]) constitute a solution for the system (CD). This solution 
is just a subcase of the previous one corresponding to the condition (1741 a) and the choice 
Q(f) = 0, P(t]) = 7]^/ 2 in ([94]). 

4. Final remarks. 

In this paper, we have obtained the infinitesimal generators which generate the Lie 
algebra of symmetries for the system ([T|) describing a planar flow of an ideal plastic 
material. The existence of generators Pi to P 5 , Di,D 2 , L, B 2 , given in equation (j4J), 
together with generator X\ of an infinite subalgebra, defined by (jSJ), were already known 
in the literature p.J. However, we have shown that the symmetry group is completed 
by the addition of generators Bi, P3, B4, B§, Bq, K, defined by (T4j) and by X 2 
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given to equation flS]), which generate an infinite-dimensional subalgebra. We have 
seen that it is possible to include the generator K in the basis of a finite-dimensional 
Lie subalgebra only if no generator B3 to B§ and Pi to P4 appears. For this reason, 
we consider separately the case of the subalgebra L which excludes K and that of 
subalgebra S including K, defined by equations (J5J) and ([9]) respectively. For each 
of these two subalgebras, we have performed a classification of the subalgebras into 
conjugation classes under the action of the symmetry group using the methods described 
in (TTJ [12] (see Section [2] and the appendix). This classification is an important tool 
in the analysis of invariant and partially invariant solutions. A classification of the 
symmetry subalgebras of (CEJ) has been performed in the past [1] for one-dimensional 
subalgebras. However, the classification presented here is more complete in the sense 
that it includes new infinitesimal generators and a classification of two-dimensional 
subalgebras, which can be used to obtain partially invariant solutions. In section |3l we 
have performed (as an example) symmetry reductions corresponding to one-dimensional 
subalgebras represented by newly found generators. For the first reduction, we have used 
the generator B\ and for the second reduction the generator K, both defined in the list 
of generators fll]). The symmetry reduction, using the invariants of B, leads to a new 
solution (see 021 p . ([22]) ) defined in terms of an arbitrary function of £ = xu + yv and 
where the velocity fields are implicitly defined. For a particular choice of the arbitrary 
function, we have traced (in Figured]) the shape of an extrusion die corresponding to this 
solution. The obtained solution has the particularity that the kinetic energy is constant 
along the flow. A similar analysis has been performed for the generator K in order to 
obtain a partially invariant solution. For this solution, the invariants given in equation 
(124]) were used to add constraints which allow us to obtain a solution more easily. These 
considerations were illustrated by finding a particular solution of system ([1]) defined by 
equations f l27|) . (128]) et fl30|) . An example of a velocity vector field and an example of 
an extrusion die have been traced respectively in figures [2] and [3j It should be noted 
that, to the vector field in Figure [2] {i.e. for the corresponding parameter values), we 
can associate a large family of extrusion dies, of which the one in Figure [3] is a particular 
choice. The contours have to be chosen in such a way that they correspond to flow lines 
to the extent that it is possible to trace two curves linking them which satisfy equation 
(T5]). These curves are the limit of the plasticity region. 

An interesting observation concerning the generator K is that if we take the 
commutator of K with the generators {P\, P2, P3, Pa}, we obtain the generators 
{B 3 , B4, B 5 , Bq}. Repeating the procedure with generators {B 3 , B4, B 5 , B 6 }, we generate 
four new generators, and so on. An interesting fact is that, at each step, the new obtained 
generators can complete the base of the subalgebra Li generated at the previous stage in 
order to form a new higher-dimensional (but still finite-dimensional) subalgebra (ensured 
by excluding K from the base). Consequently, it is always possible to enlarge a finite- 
dimensional subalgebra Li that excludes K by increasing its base with the result of 
commutators [K, Z] with Z 6 Li. This gives us a chain of finite-subalgebras of the form 
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£c£iC£ 2 C...c£,C This subject will be addressed in future works. 
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Appendix. List of two-dimensional representative subalgebras of C 



{Bi,D 2 \ 

{B u D 2 + \L + a 1 D l + a 2 P 5 } 
{B 1 ,D 2 + XD 1 } 
{B 1 ,D 2 + 3D 1 +5P 5 } 



£2,23 = {Bi + eD l + dP 5 ,L + aP b \ 
r 2 ,24 = {-Bi + e£> 1 +aP 5 ,Pi} 
C 2 ,25 = {Bi + eB 5 +aB 6 ,B 3 } 

£2,26 = {Bi + 61B5 + aiB 6 , B 3 + a 2 Pi + a 3 P 2 + e 2 P 4 } 
£2,27 ={Bi + eiB 5 + ai B 6 + a 2 P 4 + e 2 P 5 , 

B 3 + a 3 Pi + a 4 P 2 + a 4 P 4 } 
£2,28 = {Bi + eB 5 + aiP 4 , D 2 + 3P X + a 2 P 5 } 
£2,29 = {Pi + eP 5 + aiPi + a 2 P 5 , Pi} 
£2,30 = {Pi +eP 3 ,P 5 } 
£2,31 = {Pi + (P 3 + aiP 4 , P 3 + a 2 P 2 } 
£2,32 = {Pi + £iP 3 + aiP 4 + e 2 P 5 , B 3 + a 2 P 2 } 
£2,33 = {Pi + £iP 3 + aiP 4 + a 2 P 5 , P 3 + a 3 P 2 + e 2 P 4 } 
£2,34 = {Pi + eiP 3 + aiPi + a 2 P 5 , B 3 + e 2 P 4 + SP 5 } 
£2,35 = {Pi + (P 3 + aiP 4 + a 2 P 5 , P 3 + SP 5 } 
£2,3% = {Pi + (P 3 + aPi + e 2 P 5 , Pi} 

£ 2 ,37 = {Pi+EiP 3 + aiP 4 ,Pi} 

£ 2 ,38 -{Pi + eiP 4 ,Pi} 

£2,39 = {Pi + eiP 4 + e 2 P5,Pi} 

£2 4 o = {Pi + ePi, L + aiPi + a 2 P 5 } 

£2,41 = {Pi + eP 5 ,P 3 + aiP 2 + a 2 P 4 } 

^2,42 = {Pi + eP 5 ,B 3 + aP 4 + e 2 P 5 } 

>C2, 43 = {Si+eP5,Pi} 



£2,1 

£2,2 
£2,3 
£2,4 

£ 2 , 5 = {B 1 ,D 2 -D 1 + 5P 5 } 



£2,6 '- 
£2,7 - 
£2,8 = 
£2,9 = 
£2,10 
£2,11 
£2,12 
£2,13 

^2,l 4 
•^2,15 
£2,16 
£2,17 
£2,1S 
£2,19 
£2,20 
£2,21 
£2,22 



= {B 1 ,D 2 + D 1 +5P 5 } 

= {B U D 2 + 5P 5 } 

= {B 1 ,L + a 1 D 1 + a 2 P 5 } 

= {Pi, Pi +aP 5 } 

= {Pi,P 3 } 

= {B 1 ,B 3 + aP 2 + eP i } 
= {Pi,P3 + eiP 4 + e 2 P 5 } 
= {P 1 ,P 3 + eP 5 } 

= {Pi,Pi} 
= {Pi,P 5 } 

= {Pi + eL, D 1 + a x L + a 2 P 5 } 

= {B 1 +eL + aD 1 ,Ps} 

= {Pi + eL + SP 5 , Pi + aii + a 2 P 5 } 

= {Pi +eD u L + aP 5 } 

= {Pi+ePi+aiP5,P 3 + a2P2} 

= {Pi + ePi + aiP 5 ,P 3 + a 2 P 2 + eP 4 } 

= {Pi+ £ Pi,P 5 } 



Table 6. List of 2-dimensional representative subalgebras of £ that have a nonzero 
component Pi but with no P 2 component. The parameters are e = ±1 and a, A, S, e M, 
where 6 ^ 0, A > 0. 
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/* 

*- / 2,44 








£2,75 


{P 2 + Pi + eP 5 + dP 5 ,Pi} 




r 

*-'2,45 




\- L/ 2; ^1 T ui 5 j 




r 

*-2,76 — 


{D 2 + D 1 + eB 5 + 5P 5 ,P u } 




-^2,46 




IDo Bi + a Pot 

1^2? J-M ~ WJ 2j 




£2,77 — 


{P 2 + Pi + eP 6 ,P 5 } 




-^2,47 




\- L -'2j -^5 ~ WJ 4j 




£2,78 — 


{D 2 + D 1 + eB 6 ,B 5 + 6P u l 


-aP v } 


-^2,48 




IP? PA 




£2,79 


{D 2 + D 1 + eB 6 ,B 5 + SP v } 




-^2,49 




{P 2 ,P 3 } 




£2,80 — 


{D 2 + Di + eB 6 + 5P 3l B 5 + ai P 3 + a 2 PA 


-^2,50 




{P 2 ,P 5 } 




£2,81 — 


{D 2 + D 1 + eP 3 + aP i ,B 3 } 




r 

*- / 2,51 




{D 2 + XL,D 1 + aL} 




r 

^2,82 — 


W 2 + D 1 + eP 3 + aP 4 ,B 3 + 8P 2 } 


r 

- z -'2,52 




{D 2 + XL + aD 1 ,P 5 } 




f — 

-^2,83 


Wo + D^ + eP) + aPi Br\ 

\-ls 2 ' ^1 ~ CJ 3 ~ Uj± 4i ^0 J 




r 

*"2,o3 




{£> 2 -£>!,/>!} 




r 

*"2,84 — 


-f Po 4- Pi 4- eP? 4- aPi Bk 4 


SPa\ 


-^2,54 




{P 2 -Pi,P 3 + eP 5 } 




£2,85 


{^ 2 + ^ 1 + eP3 + aP4,P 5 + 


8P 5 } 


-^2,55 




{D 2 -D 1 + eB 3 ,P 1 } 




£2,86 


{ J D 2 + £) 1 + eP 3 + aP4,Pi} 




-^2,56 




{P 2 - Pi +eP 3 ,P 3 } 




£2,87 


{Dii + ^i + ePj + aPi.Ps} 




-^2,57 




{P 2 - D 1 + eP 3 + aiP 2 + a 2 P 5 , Pi} 


r 

*-2,88 — 


{D 2 + D 1 + eP 4 ,B 3 + aP 2 } 




r 

*- , 2,58 




{D 2 - D 1 +eP 1 +a 1 P 2 ,B 3 - 


^a 2 P 2 } 


r 

*-'2,89 — 


{D 2 + D 1 + eP 4 ,B 5 + aP v } 




-^2,59 




{D 2 - D 1 +eP 1 +aP 2 ,B 3 + 


SP 5 } 


£2,90 


{D 2 + D 1 +eP 4 ,B 5 + SP 5 } 




r 

*-'2,60 




{D 2 -D 1 + eP 1 +aP 2 ,B 5 } 




f 

*—2,91 — 


{D 2 + D 1 + eP 4 ,P 1 } 




-^2,61 




{D 2 - D 1 +eP 1 +aP 2 ,B 5 + 


SP4} 


£2,92 — 


{D 2 + £>i + eP 4 , P 3 } 




/* 

- i -'2,62 




{D 2 -D 1 + eP 1 +aP 2 ,P 3 } 




r — 

*"2,93 


p. + A + ^Ps} 




r 

*--2,63 




{D 2 -D 1 + eP 1 +aP 2 ,P 5 } 




r 

*-*2,94 — 


{Z) 2 + A£>i, L + aPs} 




r 

*-'2,64 




{D 2 - D 1 + eP 2 , P 3 + a 2 P 2 } 




r 

*"2,95 — 


{D 2 + XD 1 ,B 3 + aP 4 } 




-^2,65 




{D 2 - D 1 +eP 2 ,B 3 + 8P 5 } 




£2,96 = 


{D 2 + XD u B 5 + aP4} 




-^2 66 





{D 2 - D 1 + eP 2 , B 5 + aP 4 } 




>C 2 97 = 


{D 2 + XD 1 ,P 3 } 




-^2,67 




{P 2 - Pi +eP 2 ,P 3 } 




£2,98 = 


{D 2 + XD 1 ,P 5 } 




^2,68 




{P 2 - Pi + eP 2 ,P 5 } 




£2,99 = 


{D 2 + XD 1 +8P 5 ,L + aP b } 




^2,69 




{D 2 -D 1 + 5P 5 ,P 3 } 




-£2,100 = 


--{D 2 + XD 1 +aP 5 ,PA 




£2,70 




{D 2 +D 1 ,P 2 } 




^2,101 = 


-- {D 2 + XP 5 ,L + aiDi + a 2 P 5 } 


-^2,71 




{D 2 + D 1 ,B 5 + eP 5 } 




£2,102 = 


= {D 2 + XP 5 ,D 1 + aP 5 } 




-^2,72 




{P 2 + Pi+eP 5 ,Pi} 




•£2,103 = 


= {D 2 + XP 5 ,P 1 } 




-^2,73 




{D 2 + D 1 +eB 5 + SP 4 ,P 1 } 




£2,104 = 


= {D 2 + XP 5 ,P 3 } 




£2, 74 




{D 2 +D 1 + eB 5 + aP4,P u } 









Table 7. List of 2-dimensional representative subalgebras of L that have a nonzero 
component D 2 but with no component in {Pi,P 2 }. The parameters are e = ±1 and 
a, A, 5, e WL, where S ^ 0, A > 0. 



£2,105 — {Pi 
£2,106 = {Pi 



P 2 ,P 5 j 

P 2 + SL, Pi + aL} 



£2,114 = {Pi - P 2 + eiL + d"P 3 , Pi - eiP 3 + e 2 P 5 j 
£2,115 = {Pi - P2 + eiL + aP 3 + SP 4 , 

Pi -eiP 3 + e 2 P 5 } 
£2,116 - {Pi - P 2 + eiL + aiP 3 + a 2 P 4 + <5P 5 , 

Pi - £l P 3 + eP 5 } 
£2,117 = {Pi -B 2 +SD U P 5 } 
£2,118 = {Pi - P2 + 5Di + aiP 5 , L + a 2 P 5 } 

£2,119 = {Pi - P 2 + aiP 5 , L + a 2 Pi + a 3 P 5 } 

£2,120 = {Pi - P 2 + «iP 5 , Pi + a 2 P 5 } 
£2,121 = jL + aPi,P 5 } 

£2,122 = {L + aiP 5 , Pi + a 2 P 5 } 



£2,107 = {Pi - P 2 + eiL, P - eiP 3 



f e 2 P 5 } 
eP 3 } 



£2,108 = {Pi - P 2 + ei + 0D1, Pi - 
£2,109 = {Pi -B 2 + SL + 0D1, P 5 } 
£2,110 = {Pi - P 2 + + aiPi + a 2 Pi + a 2 P 2 , 

P 3 - eP 5 + a 4 P 2 + a 5 P 4 + a 6 P 5 } 
£2,111 = {Pi -B 2 + eL + aPi + aiPi + a 2 P 2 , 

P 4 -eP 6 + a 3 Pi+a 4 P 3 } 
£2,112 - {Pi -B 2 + eL + 0D1 + 5P 5 , Pi - eP 3 } 
£2,113 = {Pi ~ P2 + eiL + e 2 Pi, P5} 



Table 8. List of 2-dimensional representative subalgebras of C that have a nonzero 
component P 2 . The parameters are e = ±1 and a, A, 8, £ R, where <5 7^ 0, A > 0. 
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£2,123 


= {D U B 3 } 






£2,127 


= {D 1 ,B 3 + aP 2 + eP 4 } 


^2, 124 


= {£>i, £3 + 


eB 5 4 


ai-Pi + a 2 P 2 + 03^4} 


£2,128 


= {D 1 +aP 5 ,P 1 +eP 3 } 


•^2,125 


= {r>i,B 3 + 


eP 3 } 




£2,129 


= {D 1 +aP 5 ,P 3 } 


£2,126 


= {r>i,B 3 4 


a 1^*2 


+- eP 3 + 02^4} 





Table 9. List of 2-dimensional representative subalgebras of £ that do not have 
components in {Bi, Z? 2 , -B 2 , L} with a nonzero D\ component. The parameters are 
e = ±1 and a, A, 5, € M, where 5 ^ 0, A > 0. 



£2,130 
-£2,131 
£2,132 
-£2,133 
^2,134 
£2,135 
£2,136 
£2,137 
£2,138 
£2,139 
£2,140 
£2,141 
£2,142 
£2,143 
£2,144 
£2,145 
£2,146 
£2,147 
£2,148 
£2,149 
£2,150 
£2,151 
£2,152 
£2,153 
£2,154 
£2,155 
£2,156 
£2,157 
£2,158 
£2,159 



{B 3 ,P 1 +aP 2 \ 
{B 3 + e 1 B 5 ,P 1 +aP 2 + e 2 P 3 } 
{B 3 + e 1 B 5 ,P 1 +e 2 P 3 } 
{B 3 + e 1 B 5 ,P 2 + e 2 P 3 } 
{B 3 + eP 3 + aP 4l P 1 +5P 2 } 
{B 3 + eP 4 ,P 1 +5P 2 } 
{B 3 + a 1 P 2l P 1 } 
{B 3 + a 1 P 2 + eP 3 + a 2 P 4 ,P 1 } 
{B 3 + aP 2 + eP 4 ,P 1 } 
{B 3 + eP bl P 1 +aP 2 } 
{B 3 + 5P 2 + 01P4, Pi + e^3 + a 2 P 4 } 
{B 3 + 01P4, Pi + eP 3 + a 2 P 4 } 
{B 3 + eP 5l P 1 +eP 3 + aP 4 } 
{B 3 + SP 3 + aP 4 ,P 2 + eP 4 } 
{B 3 + SP 4 ,P 2 + eP 4 } 
{B 3 + aP 5 , P 2 + eP 4 } 
{B 3 + ai P 2l P 3 + a 2 P 4 } 
{B 3 + ai P 2 + eP 4 ,P 3 + a 2 P 4 } 
{B 3 + eP 5 , P3 + aP 4 } 
{B 3 + aP 2l P 4 } 
{B 3 + aP 2 + eP 3l P 4 } 
{B 3 + eP 5 ,P 4 } 
{B 3 + e 1 B 5 + 5P 2 ,P 1 +aP 2 } 
{B 3 + eiB 5 + ai P 2 + e 2 P 4 , Pi + a 2 P y } 
!R 3 + eB 5 + ai P 5 ,Pi+a 2 P 2 } 

+ aiP 3 + a 2 P 4 , P 1 + a 3 P 2 4 
+ a 1 P 4 ,P 1 +a 2 P y + e 2 P 3 } 
Pi + a 2 P 2 + e 2 P u } 
Pi + aP 2 + e 2 P 3 } 
+ a 1 P 3 + a 2 P 4 ,P 1 +e 2 P 3 } 



£2,163 


= 


{B 3 


£2,164 


= 


{B 3 


£2,165 


= 


{B 3 


£2,166 


= 


{B 3 


£2,167 


= 


{B 3 


£2,168 


= 


{B 3 


£2,169 


= 


{B 3 


£2,170 


= 


{B 3 


£2,171 


= 


{B 3 


£2,172 


= 


{B 3 


£2,173 


= 


{B 3 


£2,174 


= 


{B 3 


£2,175 


= 


{B 3 


£2,176 


= 


{B 3 


£2,177 


= 


{B 3 


£2,178 


= 


{B 3 


£2,179 


= 


{B 3 


£2,180 


= 


{B 3 


£2,181 


= 


{B 3 


£2,182 


= 


{B 3 


£2,183 


= 


{B 3 


£2,184 


= 


{B 3 


£2,185 


= 


{B 3 


£2,186 


= 


{B 3 


£2,187 




f 7~> 

{B 3 


£2,188 




{B 3 


£2,189 




{B 3 


£2,190 




{B 3 


£2,191 




{B 3 


£2,192 




{B 3 


£2,193 




{B 3 



+ eB 5 + 6P 5 + aP 4 ,P 2 \ 

+ eB 5 + aP 4 ,P 2 } 

+ eB 5 + aP 4 , P 2 } 

+ eB 5 + 5P 3 + aP 4 , P 2 + e 2 P 3 } 

+ e 1 B 5 + SP 4 ,P 2 +e 2 P 3 } 

+ e 1 B 5 + SP 5 ,P 2 +e 2 P 3 } 



{B 3 + e 1 B 5 + SP 2 
{B 3 + e 1 B b + 5P 3 
{B 3 + e!B 5 + 5P 4 
{B 3 + e 1 B b + 8P b 
{B 3 + e 1 B 5 + SP 2 



e 2 P 3 } 



£2,160 = {B 3 + eiB 5 + 5P 3 + aP v ,P 1 + e 2 P u } 

£ 2 .i6i - {B 3 + eiB 5 + 8P 4 ,P 1 + e 2 P 3 } 
£2,162 = {B 3 + exPs + SP h ,P 1 + e 2 P u } 



eB 5 + aP 4 ,P 3 } 
eB 5 + SP 5 ,P 3 } 
aiPx +a 2 P 2 ,B 3 + a 3 Pi 



eP 3 + a 4 P 4 } 



£2,194 
£2,195 



+ aiPi + a 2 P 2 + eP 3 + a 3 P 4 , B 4 + a 4 P x + a 5 P 3 + a 6 P 4 } 
+ a 1 P 1 + a 2 P 2 + eP 4 , B 4 + + a 4 P 3 + a 5 P 5 } 
+ a 1 P 1 + a 2 P 2 , B 4 + eB 6 + a 3 Pi + a 4 P 2 + a 5 P 3 } 
+ a 1 P 1 + a 2 P 2l B 5 + a 3 B 6 + a 4 P 4 } 
+ 01P1 + o 2 P 2 , B 5 + a 3 B 6 + eP x + a 4 P 2 + a 5 P 5 } 
+ aiP x + a 2 P 2 ,B 5 + a 3 B 6 + eP 2 + a 4 P 4 } 

+ aiPi + a 2 P 2 + eP 3 + a 3 P v , B 5 + a 4 B 6 + a 5 P 1 + a 6 P 2 + a 7 P 4 } 
+ a 1 P 1 + a 2 P 2 + eP 4 , B 5 + a 3 B 6 + a 4 Pi + a 5 P 2 + a 6 P 4 } 
+ 01P1 + a 2 P 2 + eP 3 + a 3 P 4 , B 6 + a 4 P x + a b P 2 + a 6 P 3 } 
+ a 1 P 1 + a 2 P 2 + eP 4 , B 6 + a 3 P^ + a 4 P 2 + a 5 P 3 } 
+ a 1 P 1 + a 2 P 2 , B 6 + ePi + a 3 P 2 + a 4 P 3 } 
{B 3 + a 1 P 1 + a 2 P 2 , B 6 + eP 2 + a 4 P 3 } 
+ aiPi + a 2 P 2 ,B 6 + a 3 P 3 } 
+ a 1 B 4 + a 2 Pi + a 2 P 2 + eP 3 + a 3 P 4 , B 5 + a 4 P x 
+ a x B 4 + a 2 P x + a 3 P 2 + eP 4 , B 5 + a 4 P 1 + a 4 P 2 
+ aiB 4 + a 2 Pi + a 2 P 2 , B 5 + eP 4 + aP 2 + aP 4 } 
+ a x B 4 + a 2 P x + a 2 P 2 , B 5 + eP 4 + a 3 P 4 } 
+ a x B 4 + a 2 P x + a 3 P 2 , B 5 + a 4 P 4 } 
+ eB 5 + aiPi + a 2 P 2 + a 3 P 3 + a 4 P 4 , 

B 4 + a 5 B 5 + a 6 B & + a 7 P x + a$P 3 + a 9 P 4 } 
{B 3 + eB 5 + 01P1 + a 2 P 2 + a 3 P 3 + a 4 P 4 , B 5 + a 5 B 6 + a 6 Pi + a 7 P 2 4 
{B 3 + eB 5 + aiP + a 2 P2 + a 3 P 3 + a 4 P 4 , B 6 + a 5 P 1 + a 6 P 2 + a 7 P 3 } 



+ C15P2 + a e P 4 } 
+ a 5 P 4 } 



a 8 P 4 } 



Table 10. List of 2-dimensional representative subalgebras of £ that do not have components in {B x , D 2l B 2 , L, D±} with a nonzero P 3 
component. The parameters are e = ±1 and a, A, 5, £ R, where S ^ 0, A > 0. 



Group analysis of an ideal plasticity model 



32 



£2,196 — {-85, Bq\ 
£2,197 = {B 5 ,P 3 + aP 4 } 
£2,198 = {B5, P4} 

£2,199 = {B 5 + eP 2 + aiPi, Pi + a 2 P 2 } 

£2,200 = {B 5 + 01P4, Pi + C12P2} 

£2,201 = {B 5 + eP 5 ,P 1 +aP 2 } 

£2,202 = {B 5 + eiP 2 + oiP 4) Pi + a 2 P 2 + e 2 P 3 } 

£2,203 = {B 5 + 01P4, Pi + «2P 2 + eP 3 } 

£2,204 - {P 5 + eiP 5 , Pi + aP 2 + e 2 P 3 } 

£2,205 = {P 5 + ePi+aP 4 ,P 2 } 

£2,206 = {B 5 + aP 4 ,P 2 } 

£2,207 = {B 5 + eP 5 ,P 2 } 

£2,208 = {B 5 + SJ\ + «P 4 , P 2 + eP 4 } 



{P 5 +aP4,P2+eP 4 } 

{P 5 +'5P5,P2 + eP 4 } 

{P 5 + ePi + aip2 + a 2 P 4 , P3 + a 3 P v } 

{B 5 + eP 2 + aiPi,P 3 + a 2 Pi} 

{B 5 + 6P 4 ,P 3 + a 2 P 4 } 

{B 5 + eP 5 ,P 3 + aP 4 } 

{B 5 + eP 1 +aP 2 ,P 4 } 

{B 5 +eP 2 ,P 4 } 

{B 5 +eP 5 ,P 4 } 

- {B 5 + a x P 3 + a 2 P 4 , 

B e + ePi + a 3 Pi + a 3 P 2 + a 4 P 3 } 
{B 5 + ai P 3 + a 2 P 4 , B 6 + eP 2 + a 3 P 3 } 

- {B 5 + ePi + aiP 2 + a 2 P 3 + a 3 P 4 , 

B e + a 4 Pi + a 5 P 2 + a 6 P 3 } 
{P 5 + eP 2 + aiP 3 + a 2 P 4 , P 6 + a 3 Pi + a 4 P 2 | 



£2,209 : 

£ 2 ,210 : 
£2,211 : 
£2,21 2 : 
£ 2 , 213 : 
£2,21 4 : 
£2,215 : 
£2,216 : 
£2,217 : 
£ 2 , 2 18 

£2,219 : 

£ 2 , 22 o 



-2.221 = 



Table 11. List of 2-dimensional representative subalgebras of C that do not 
have components in {Pi, D 2 , B 2 , L, D\, B 3 , B 4 } with a nonzero P5 component. The 
parameters are e = ±1 and a, A, S, <G R, where S 7^ 0, A > 0. 



£ 2 . 222 — {Pi, P5} 

£2.223 = {P +aP 2 ,P 5 } 

£ 2> 2 24 = {Pi + aP 2 + eiP 5 , P 2 + e 2 P 5 } 

£2.225 = {Pi +eP 3 ,P 5 } 

£2,226 = {Pi + eiPs + e 2 P 5 , P 2 + aiP 3 + a 2 P 4 + SP 5 } 
£2,227 = {Pi + eiP3 + e 2 P 5 , Pv + 6P5} 
£2.228 = {Pi+eiP 5 ,P2+e 2 P 4 } 
£2.229 = {Pi + eP 5 ,P 2 +6P 5 } 
£2,230 = {Pi +eP 5 ,P 3 +5P 4 } 

£2.231 = {P +£iP 5 ,P3 + e 2 P 5 } 



£2,232 = {P +eP 5 ,P 4 } 

£2,233 = {P2, P5} 

£2,234 = {P 2 + aiPs + a 2 P 4 , P 5 } 
£2,235 = {P 2 + eP5,P3 + e 2 P5} 

£2,236 = {P2 + eiP5,P3 + e 2 P 5 } 

£2,237 = {P3, P5} 

£2,238 = {P 3 + aP 4 ,P3} 

£2,239 = {P 3 + eP5,P 4 + 5P 5 } 
£2,240 = {P4, P5} 



Table 12. List of 2-dimensional representative subalgebras of C corresponding to 
translations with the P5 component. The parameters are e = ±1 and a, A, 5, € M, 
where 5 ^ 0, A > 0. 



£ 2 , 2 41 




£ 2 , 246 


= {Pi+eP 3 ,P 2 4 


oiP 3 + a 2 P 4 | 


£ 2 , 2 4 2 


= {P,P 2 +eP 4 } 


£2,247 


= {P 1 +eP,,P 3 4 


aP 4 } 


£ 2 , 2 43 


= {Pi,P 3 + aP 4 } 


£ 2 , 248 


= {Pi+eP 3 ,P 4 } 




£ 2 , 2 44 


= {Pi,P 4 } 


£ 2 ,249 


= {P 2 ,P 3 } 




£ 2 , 2 45 


= {P 1 +aP 2 ,P 3 } 


£2,250 


= {P 3 ,P 4 } 





Table 13. List of 2-dimensional representative subalgebras of C corresponding to 
translations without the P 5 component. The parameters are e = ±1 and a,X,S,£ R, 
where S ^ 0, A > 0. 



